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Abstract 

CN ■ 

We consider the (2 + 1) flavor Polyakov Quark Meson Model and study the fluctuations (correlations) of con- 
served charges upto sixth (fourth) order. Comparison is made with lattice data wherever available and overall 
i-C ' good qualitative agreement is found, more so for the case of the normalised susceptibilities. The model predictions 

for the ratio of susceptibilities go to that of an ideal gas of hadrons as in Hadron Resonance Gas Model at low 
temperatures while at high temperature the values are close to that of an ideal gas of massless quarks. Our study 
provides a strong basis for the use of PQM as an effective model to understand the topology of the QCD phase 
diagram. 

^ ! PACS numbers: 05.40.-a, 12.38Aw, 12.38Mh, 12.39.-x 

>: 

i/-) 1 Introduction 

CO; 

CO | The QCD medium is known to undergo phase transitions [1,2]. The QCD Lagrangian for Nt flavors of massless 
■ quarks has a SUl (Nf) x SUr (Nf) global symmetry which spontaneously breaks into SUy (Nf) in the low energy 
f^) \ hadronic vacuum by the formation of chiral condensate together with Nj — 1 massless Goldstone bosons. In the 
CNI ■ opposite limit of infinite quark mass, QCD becomes a pure SU (3) gauge theory. The low energy vacuum that 
' m 1 1 possesses a center symmetry Z (3) under the color gauge group (confined phase) gets spontaneously broken in the 
^ ■ high temperature/baryon density regime (deconfined phase). The study of the various phases of strongly interacting 
matter has been the subject of intense theoretical research for some time now. With data coming in from RHIC and 
LHC, experiments too are not far behind and will provide us a unique opportunity to improve our understanding of 
cd | the thermodynamics of strongly interacting matter. 

At these temperatures and densities QCD is a strongly coupled theory and hence perturbative methods to study 
the details of this phase transition fail. First principle Lattice QCD (LQCD) Monte Carlo simulations give us im- 
portant insights into various aspects of the phase transition. However LQCD suffers from the sign problem at non 
zero baryon density [3,4] and although several methods have been developed [3,4] to bypass the sign problem at 
small baryon chemical potential, a satisfactory solution for all values still eludes us. Another approach is to study 
various phenomenological models whose phase diagrams possess the essential features of QCD. These models serve to 
complement LQCD computations and also provide us with an intuitive and physical understanding about the behavior 
of the phases of strongly coupled matter in regions that are both accessible and inaccessible to LQCD. 

QCD has three conserved charges: baryon number B, strangeness S and electric charge Q. Thus the thermodynamic 
state of the bulk medium can be specified by four quantities: temperature T and the chemical potentials corresponding 
to the conserved charges namely /is, fiQ and us. Effective QCD-like models show first order hadron-QGP phase 
transition for large fig. This line of first order phase transition is expected to end at a critical end point (CEP) 
at finite \xb since general symmetry arguments [5] as well as LQCD [6-9] indicate a smooth crossover at \ib = 0. 
Divergent fluctuations and correlations are characteristics of critical behavior. In QCD, fluctuations and correlations 
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of the conserved charges are expected to provide telltale signatures of the CEP [10] and bring forth the structure of 
the QCD phase diagram. These fluctuations can also be extracted experimentally through event by event analysis 
making them important observables in understanding the nature of strongly interacting matter [11,12]. 

There have been numerous studies of the QCD correlators, both in lattice as well as models. LQCD with 2 
flavors found the fluctuations to rise steeply [13, 14] while the higher order susceptibilities exhibited non-monotonic 
behaviour [15, 16] in the crossover region. Similar studies have been performed in the 2 flavor Polyakov Nambu 
Jona Lasinio (PNJL) model [17-21], 2 flavor Quark Meson (QM) model [22], 2 flavor Polyakov Quark Meson (PQM) 
model [23] as well as in Hard Thermal Loop approximation [24-26] . Susceptibilities were also computed in an improved 
version of the 2 flavor PQM model with the addition of the frequently neglected fcrmion vacuum fluctuations [27,28]. 
For (2 + 1) flavors, there are LQCD measurements [29-32] as well as model studies in PNJL [33-37] and PQM [38,39]. 
Recently, higher order generalised susceptibilities and their sign structure in the phase diagram was studied in the 
(2 + 1) PQM model [40]. 

In this paper, we study the diagonal and off-diagonal susceptibilities of the conserved charges B, Q and S in the 
(2 + 1) PQM model at zero chemical potentials with the inclusion of the contribution from the fermionic vacuum 
fluctuations. The significance of the vacuum term was recently pointed out in the case of 2 flavours [27] where it was 
shown that the order of the phase transition in the massless limit changes from first order to crossover on adding the 
vacuum term. In our earlier study [41], we have shown the significance of this vacuum term on the thermodynamics in 
the case of (2 + 1) flavors. The phase transition region became smoother which led to better agreement with lattice. 
Here we report our investigation on the susceptibilities of this model and how they compare with results from the 
HotQCD [31] and Wuppertal-Budapest (WB) [32] LQCD groups. Our paper is organised as follows: In Sec. 2, we 
discuss the details of the PQM model and its parameters. In Sec. 3, we present our results and compare with those 
obtained in LQCD. Finally, in Sec. 4 we summarise and conclude. 



2 FORMALISM 



2.1 Model 



We had formulated the (2 + 1) flavor PQM with the inclusion of the vaccum term in [41]. Here we work with the 
same model. The relevant thermodynamic potential at a temperature T and chemical potentials ^s, fiQ and /xg in 
the mean field approximation can be written as [41] 

(T,HB,HQ,HS) = U M (<T u ,<T d ,(T s ) +Wp ly-VM ($,$,T) + ilgg (T, fi B , fiQ , [i S , , <?d, &s , $) (1) 

Here 07 denotes a chiral condensate of the quark with flavor / (up(u), down(d) and strange(s)). Um {&u,&d, o"s) is the 
contribution from the mesonic potential given by 

U I \ ^ f 1 I 2 , 2\ , 2\ h u h d C 

UM\Vu,Vd,Vs) = — I 2 { a u + a d) + a s J - ^" cr « - -^°d - h s a s - -^-j=a u a d (T s 

+ — (°u + a d) a s H ' u +(J d+ 4cr s j + Y autTd ( 2 > 

The Polyakov loop potential Up \y-VM ($, with the Jacobian Van der Monde term is given by [42] 

^Poly-VM W Poly ($,!>, T) r - , 3 _ 3 . . 

— — \ '- = y -±~i '- - k log 1 - 6$$ + 4 ($ J + $ J ) - 3 ($$) 



rp4 rp4 

Here Up Q \ y is the Ginzburg-Landau type potential given by [43 



(3) 



T 4 2 6 K 'A 

with 

To , fT \ 2 , /V 3 



b 2 (T) = ao + ai ^+a 2 ( ^1 +a 3 (^) (5) 
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where 



a = 6.75 , ai = -1.95 , 6 3 = 0.75 
a 2 = 2.625 , a 3 = -7.44 , b 4 = 7.5 

The quark/ antiquark contribution is given by 

ttqq (T,fl B ,HQ,Hs,Vu,Vd,<Ts,®, $) = Wq q (<T U , <T d , <J S ) + fi)£ (T ', fl B , MQ, MS , <T U , <?d, O s , $) 



(6) 



is the quark/ antiquark contribution due to thermal fluctuations 



ttfq (T,(j,B,HQ,Hs,<ru,<rd,<rs,®,®) 



-2T 



£ / 



f= u,d,s 



-2T J2 

f=u,d, 



d 3 p 

(2tt) 3 

d 3 p 
(27)3 



In 



In 



1 + 3<f> e - £ / +/T + 3$ e - 2B / +/T - 
1 + 3i>e~ £ / /T + 3$e~ 2B 7 /T 



e -3J5+/T 
. e SEJ/T 



(7) 
(8) 



where 

M/ denote the quark chemical potentials fi u , pd and Ms- They are related to mbj fiQ and by the following 
transformations 



Mm = 


Ms 
3 


, 2 ^Q 
3 


Md = 


Ms 
3 


MQ 
3 


Ms = 


Ms 
3 


3" 



In (8), is the single particle energy of a quark/antiquark. 

Ef = \Jp 2 + mj 2 
and m/ is the mass of a quark/antiquark with flavor / given by 



(9) 
(10) 



771/ 



-07, / eu,d 



V2 



07, / e s 



(11) 



fi| 9 is the vacuum term 



(CT u ,<7d,CT s ) 



2iV c 2 / 

f—u,d,s 



d 3 p 



Nr. 



m f 



A 



(12) 



f—u^d.s 

where A is the regularisation scale. In [41], we have shown that at the mean field level, f2 and hence all the thermo- 
dynamic quantities are independent of the choice of A. The model parameters: m 2 , Ai, A2, c, g, h u , hd and h s are 
determined by establishing the vacuum properties. The parameters of ^(T), 63, 64, To and k are obtained by fitting to 
the lattice data [41,43]. Two sets of model parameters were obtained in our earlier work [41], namely ModclHotQCD 
and ModelWB. Here we work with the same model parameters. In Table 1, we list the values of all the parameters 
that have been used in this work. 

The mean field values of the condensates a u , Od and <r s and the Polyakov loop variables $ and <l are determined 
by numerically solving the following set of equations 

dfl _ dfl dtt 
du u d<7 d da s 



dn on _ 



(13) 
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Model 


m 2 [MeV 2 ] 


Ai 


A 2 (A) 


c 


h u = h d [McV 3 ] 


h s [McV 3 ] 


9 


T 


K 


ModelHotQCD 


80647.587 


-8.165 


138.45 


4801.95 


1.785 x 10 b 


3.805 x 10 Y 


6.5 


210 


0.1 


ModelWB 


80647.587 


-8.165 


138.45 


4801.95 


1.785 x 10 B 


3.805 x 10 Y 


6.5 


270 


0.2 



Table 1: The parameter sets obtained with A = 200 MeV. 



X 


B 


Q 


S 


xi 


1/3 


2/3 


1 


Xa 


2/9^ 2 


4/3^ 2 


6/tt 2 


XY 


BQ 


BS 


QS 


xtt 





-1/3 


1/3 


xtt 


4/9tT 2 


2/37T 2 


2/37T 2 


XY 

A31 





-2/97T 2 


2/9tt^ 


XY 
A13 


4/97T 2 


-2/tt' z 


2/tt' z 


XYZ 


BQS 


QSB 


SBQ 


XYZ 
A211 


2/97T 2 


-2/9tt 2 


-2/3tt 2 



Table 2: The SB values of the susceptibilities and correlations. 



2.2 Susceptibilities 

The pressure P is given by 

P (T, fj, B , (j,q,hs) = -fi (T, [is-, HQ, Hs) (14) 
In order to find the generalised susceptibilities of the conserved charges, one has to take appropriate derivatives of P 

BQS = d^+HPjT*) 

Xl]k d^B/TYd^Q/TY d([i s /T) k 

In this work, we perform all our computations at zero chemical potentials. The baryon number, electric charge, and 
strangeness densities vanish at zero chemical potentials while the higher order derivatives with (i + j + k) even are 
nonzero. In this work, we compute xfjk S u Pto (i + j + k) = 6 order. The derivatives in (15) have been computed 
using the algorithmic differentiation techniques available in ADOL-C [44] which allow us to compute higher order 
derivatives efficiently and without additional truncation errors. 



3 Result 



In the low T regime, the condensates have large values which impart large mass to the relevant degrees of freedom. 
Thus all susceptibilities are small in the low T regime. In the crossover regime, the condensates start to melt and thus 
the relevant degrees of freedom become lighter. Generically, this results in a higher value of the susceptibilities as one 
approaches T c . These quantities have been extensively computed on the lattice [13-16,29-32]. In the high temperature 
limit when T is the only relevant scale, one expects the system to behave like a Stefan-Boltzmann (SB) gas of 3 flavors 
of massless quarks ( we don't consider the gluons as they do not carry any B, Q and S quantum numbers and hence 
don't contribute to the susceptibilities in the high temperature limit ) whose pressure is given by 



pSB 

rp4 



f=u,d, 



60 



+ 



2 V T J 



1 

4V 2 



T J 



(16) 



The corresponding SB limit values for quadratic and quartic fluctuations of conserved charges as well as their 
correlations are given in Table 2. On the other hand, at the low T limit, an ideal gas of hadrons as in hadron 
resonance gas models (HRGM) is expected to describe the thermodynamics and susceptibilities. Latest LQCD data 
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from WB [32] exhibit the above expectations quite clearly. Since the PQM model has quark degrees of freedom, in 
the the high T limit, one expect that the susceptibilities will go to that of an ideal SB gas of quarks. On the contrary, 
in the low T side, since the PQM model does not have baryons (there are only 'three quark' states which mimic 
baryonic degrees of freedom in the low T limit in the model [48]), one should not expect the model predictions of the 
susceptibilities to match with those of HRGM. However, in case of ratios of susceptibilities, we argue below that that 
in the low T limit the PQM results match with that of HRGM [48]. 

Let us first look at the susceptibility ratios for low temperatures in PQM. At low T, <p and (j> are nearly and 
hence from (7) we find that the excitations of one and two quark states are suppressed. Thus, we may approximately 
write the thermal quark-antiquark contribution to the pressure as [48] 



T 4 27T 2 \ T ) \ T I V T 

which is same as that of a non- interacting gas of paticles and anti-particles with baryon number 1 and -1 respectively 
and mass 3m 9 . Here d q is the effective degeneracy factor which depends on the particular susceptibility one is looking 
at: d q is 2 when the contribution from the non strange quark sector is to be considered while d q is 1 where the strange 
quark sector is important. Thus we find in the low T limit there is a nice factorization and one could write (17) 

Pqq i T ^) / TTlq \ / 3/U, 



rp4 



/(f) cosh m\ (18) 



Now we shall see that in the low T limit, a similar factorization also occurs in the HRGM framework whose pressure 
is given by that of an ideal gas of hadrons 

P HRG /r 4 --^ Yl ^Zi(T,V^ B ,IMQ,tM S ) (19) 

where 



z£ hadrons 



InZ^^j dpp 2 In (l + az t e-^ T ) (20) 

is the contribution from the ith hadron with a — =pl depending on whether it is meson or baryon respectively, 
Ei = \/ p 2 + m 2 where m; is the mass of the ith hadron with degeneracy factor di and fugacity 

z t = exp {{Bins + QiVQ + Shis) /T) (21) 

The sum in (19) is usually over all hadrons from the particle data book with masses rrii < 2.5 GeV. However, recently 
the significance of the yet undiscovered heavier hadrons on the heavy ion phenomenology and freezeout conditions 
were pointed out in [45,46]. Expanding the fugacity term one can write (20) in the following series form 

In Z t = —d t (f ) r 2 ^ (Imi/T) cosh (/ (B if x B + Q it i Q + S if i s ) /T) (22) 

Thus (19) becomes 

di 



„HR.G 



l TA = E ^ (y) E (Imi/T) cosh (/ (Bi/iB + Q l ^ Q + Sins) /T) (23) 



i=i 



Note that since —a = 1 for a meson, all the terms of the series in (23) are positive while in the case of a baryon, it is 
an alternate series. For sufficiently small temperature compared to the mass of the hadron one can approximate the 
series in / by its leading term which is the Boltzmann approximation. Also the dominant contribution will come from 
the lightest relevant hadron hi, 

p HRG /T 4 M 4 (^) 2 K 2 (f ) COSh {{Bifl B + QlfJ-Q + S&s) IT) (24) 
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Thus, once again we find a factorization in the low T limit in HRGM similar to (18) for the case of PQM. This 
factorization leads to the fact that in the low T limit ratios of susceptibilities are independent of the details of the 
mass spectrum both in the case of HRGM as well as PQM. Thus as long as the relevant degrees of freedom in both 
PQM and HRGM have the same quantum numbers, their ratios of susceptibilities should match. 

It is easy to check from (24) that ratios like xf Ixf become simply X % t ~ 3 where X e B,Q, S and Xi are the B, Q, S 
quantum numbers of hi. Since the lightest hadrons in any sector are singly charged, X] -1 is just unity and independent 
of the details of the hadron spectrum. However, in case of S and Q sectors, the contribution from multiply strange or 
charged hadrons enhances the higher order fluctuations compared to that of the second order and leads to deviation 
from unity for slightly higher temperatures. The light mass of the pions result in the breakdown of the Boltzmann 
approximation at slightly higher temperatures which in turn affects the ratios in the Q sector. As seen from (23), in 
the case of pion which is a meson, the quantum corrections are all positive and higher the order of the susceptibility 
larger is the correction. This enhances ratios likes X?/x? an d Xe /X2 '■ 

We will now present the generalised susceptibilities and correlations computed in PQM at zero chemical potentials. 
We will also present the LQCD data from HotQCD and WB groups for comparison. We find overall good qualitative 
agreement between model predictions and LQCD. Particularly in the case of ratios, the PQM predictions interpolate 
between the low temperature HRGM values and the high temperature SB values. 

3.1 (i + j + k) = 2 

There are six second order susceptibilities: three diagonal susceptibilties X2 ■> X2 \ X2 ancl three off-diagonal ones, Xi? 1 
X11 and X11 ■ All the second order susceptibilities, which measure fluctuation of the conserved charges are small in 
the low T regime and rise monotonically across the crossover temperature before saturating at high temperatures. 
This generic behaviour is seen both in the lattice data as well as in model calculations. We have plotted all these 
quantities in Fig 1. The HotQCD data for X2 > X2 an d xf indicate a steep rise across the transition region attaining 
the corresponding SB limits by T <~ 1.7 T c . On the other hand, the WB data rises gently across T c attaining 90% 
of the SB limit by T ~ 2.4 T c . For X2 > ModclWB shows very good agreement with the WB data. On the other 
hand, the rise in ModelHotQCD though faster than ModclWB, is not strong enough and fails to catch up with the 
HotQCD data. The ModelHotQCD reaches 93% of the SB limit by T - 2.4 T c . Similar high T behaviour is observed 
for X2 ano - X2 • At low temperature, X2 from WB increases appreciably as iV T is increased and the continuum limit is 
approached. This behaviour can be attributed to the fact that the dominant contribution to X2 comes from the pion. 
For the same reason, PQM at the mean field level fails to reproduce the WB data. Thus, it is necessary to incorporate 
pionic fluctuations in our present study for a better estimate of X2 (f° r similar studies in the context of the PNJL 
model, see [47]). 

The X11 correlation exhibits a peak in the transition region. As already mentioned, in the low T side because 
of the high mass of the relevant degrees of freedom, it goes to zero, while in the SB limit it is clear from (9) that 
the contribution to X11 from the u quark sector is neatly cancelled by the sum of the contributions from the d and 
s sectors and hence goes to zero again. ModelHotQCD agrees reasonably well with the HotQCD data. In the peak 
region, it is about 40% larger than ModelWB. In case of the other correlators xfi an d X11 1 the models compare in 
the same way to the lattice data as the diagonal susceptibilities. 

We shall now focus on the ratio of the above susceptibilities. We have plotted these ratios in Fig. 2. The normal- 
ized correlations appear to show a better agreement with lattice data as compared to the unnormalized correlations 
themselves of Fig 1. In these plots, the ModelHotQCD and ModelWB lie almost together. Xi? IX2 shows a peak in 
the transition region and both in the low and high T regimes, this ratio goes to zero. In the high T limit, X11 itself 
goes to zero and thus the ratio goes to zero. In the low T limit, the dominant contribution to Xi? is from the effective 
'3 quark state' from the non strange sector while the dominant contribution to xi comes from the pion. The three 
quark state being much heavier than the pion, Xi? approaches zero much faster than X2 ano - thus the ratio goes to 
zero in the low T limit. In the case of X11 /xfi there is no peak like structure and the ratio decreases monotonically. 
In this case the high T behaviour is same as Xi? IX2 ancl it goes to zero as Xi? itself goes to zero. In the low T 
limit, both X\? an d X2 S e t contribution from the '3 quark state' from the non strange sector. It is straightforward 
to check from (9) that because of the relative sign between \ib and /j,q in the expression of fid, Xi? IX2 approaches 
1/2. While — xfiVxf [12] shows a steep rise and reaches the SB limit by T/T c ~ 1.3, — XuVxif r i ses gently before 
meeting the SB limit at T/T c ~ 1.8. In the low T limit, we expect the ratios to go to zero: since the three quark states 
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0.35 




Figure 1: Plots of all the second order susceptibilities in ModelHotQCD (in red) and ModelWB (in blue). The lattice 
data in black and red are obtained using p4 action with N T = 4 and 6 respectively by the HotQCD collaboration [31]. 
The blue lattice data are continuum estimate from WB collaboration using stout action [32] . In case of x*2 > we have 
also plotted the WB lattice data for N T — 8(brown), 12(orange) and 16(green). In all the curves the high T SB limit 
is indicated by a thick black line. 
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0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

T/T c T/T c 

0.35 ^ 1 1 , 1 , 1 1 1 1 1 1 1 1 , 1 , 1 1 , ^ f - ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

T/T c T/T c 



Figure 2: Plots of second order correlations normalised to the second order susceptibilities in ModelHotQCD (in red) 
and ModclWB (in blue). The lattice data in black and red are obtained using p4 action with N T = 4 and 6 respectively 
by the HotQCD collaboration [31]. The blue lattice data are continuum estimate from WB collaboration using stout 
action [32] . In all the curves the high T SB limit is indicated by a thick black line while the dashed black line at low 
T shows the HRGM value. 
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from the strange sector (which dominantly contribute to \ii m PQM) are much heavier compared to the three quark 
states from the non strange sector and the kaons which dominantly contribute to X2 an d xf respectively. As one can 
verify from Table 2, Xii IX2 an d X11 /xf approach 1/2 and 1/3 in the high T limit by 1.8T C and 1.1T C respectively. 
Their low T behaviour is quite distinct. Since the kaons which contribute dominantly to x?i are much heavier than 
the pions that has a leading contribution to Xsn xfi /X2 S oes to zero. On the other hand, kaons alone contribute 
dominantly to both xfi and xfi an d the ratio approaches 1/2 as only two out of the four kaons are charged. X\i Ixi 
also exhibits interesting feature in the transition regime, a plateau structure near T c . Similar observation has been 
found even on the lattice [31] as well as in the case of PNJL [35,37] model. In [37] this has been attributed to the 
shift in the role of the dominant degrees of freedom from the hadrons to quark quasiparticles just above T c . 

At low temperatures, as argued earlier the ratios obtained in PQM are found to compare well with that of HRGM 
which are also shown for comparison. The correlator Xi? project only on the charged baryon sector. Thus within 
HRGM at low temperatures, the ratio xfi IX2 approach \ as X2 gets contribution from both protons and neutrons 
while xfi gets contribution only from protons. For similar reasons, X11 /xf a ^ so approach i. We note that the 
HRGM plot of X11 /xf m Fig. 2 shows a value slightly greater than 0.5 which is due to the small mass difference 
between the different isospin states of kaon (which is the lightest relevant hadron in this case) . X11 /xf g° to zevo at 
low T because the lightest baryons do not carry strangeness. Similarly, Xi? IX2 ■• XuVxf > X11 /X2 an d X11 /X2 S° 
to zero too. 

Thus we see that all the second order fluctuations increase monotonically across T c . In the case of the correlations, 
the story is more varied: while xbq exhibits a peak in the transition regime, xbs and XQS rise monotonically across 
T c . When we consider the ratios, we have an even more varied features. xfi IX2 an d ~Xii IxE show a monotonic 
beahviour, xffVxf and X\\ IX2 exhibit sharp changes in the transition regime, Xi? /X2 shows a sharp peak and 
X11/X2 exhibit a plateau in the transition regime. Thus quantities like xfi Q , Xn Q /x?, xfiVxf . Xii/xi and X11/X2 
are well suited to probe the QCD phase transition because of their distinct features in the transition regime. In fact, 
it was already pointed out in [12] that the ratio XuVxf has different values in the hadron and QGP phases and thus 
is a very sensitive tool to probe the nature of the strongly interacting matter. 

3.2 (i+j + k)=4 

The higher order susceptibilities are more sensitive to fluctuations in a system. We show in the left column of Fig 
3 plots of fourth order diagonal susceptibilities x%> Xq and x%- Unlike the second order diagonal susceptibilities, in 
this case ModelHotQCD and ModelWB almost lie on top of each for the entire temperature range. A notable feature 
for all the three susceptibilities is that they all peak around the cross-over temeperature which is also evident from 
LQCD. However, when compared to lattice the models show a smaller peak. But here we should note that the data 
for N T = 6 are consistently smaller than the N T — 4 data and so one would expect that the continuum estimates of 
LQCD will further come down which will be in better agreement with the PQM predictions. Xs shows a very gentle 
peak. This is a common feature of effective models (also seen in PNJL [35] models), where the strange quark mass 
reduces at a slower rate as compared to the lighter quarks in the transition region. In the high T region, both model 
as well as LQCD seem to attain the SB value in the case of x% an( ^ Xq while for x| the models arc 11% higher than 
the SB limit, whereas lattice is 50% higher at T/T c — 2.4. 

In the right column of Fig 3, we have shown plots of the ratio of the fourth order susceptibilities to those of the 
corresponding second order ones. As in the case of second order correlations, it is worthwhile to note that the models 
are in better agreement with lattice studies when the ratios of the susceptibilities are compared. In all the cases, 
ModelHotQCD predictions are slightly lower as compared to that of ModelWB. As argued earlier, at low T, ratios 
of the kind xf Ixf where X € B,Q,S approach unity both in PQM as well as in HRGM. In the crossover region, 
all the ratios exhibit a peak, the most prominent being that of Xa I X? which starts deviating from unity only around 
T <~ 0.8T C making it a good candidate to probe the crossover region [48]. In the high T limit, the model predictions 
saturate around the SB values. 

For completeness in Fig 4, we have presented the model computation for all the fourth order correlations. There 
are total 16 of them. LQCD data for these correlations are still not available. All the fourth order correlations stay 
at zero for temperatures upto ~ 0.7T C , while in the high T limit they all approach within 10% of the SB values. All 
the correlations in the BQ sector exhibit a sharp peak in the transition regime, while the introduction of strangeness 
broadens the peak in the other cases due to the slow melting of the strange condensate. For some cases like X^' X31 \ 
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Figure 3: On the left are plots of fourth order susceptibilities and on the right are plots of these susceptibilities 
normalized by their second order susceptibilities in ModelHotQCD (in red) and ModelWB (in blue). The lattice data 
in black and red are obtained using p4 action with N T — 4 and 6 respectively by the HotQCD collaboration [31]. In 
all the curves the high T SB limit is indicated by a thick black line while the dashed black line at low T shows the 
HRGM value. 
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Xzfi ~xii an d ~Xi2i there is a hint of a broad double peak in ModelWB, although in case of ModelHotQCD the 
second peak looks sufficiently suppressed. Similar double peak structures have also been seen in the case of (2 + 1) 
PNJL model [37]. 

3.3 (i+j + k) = 6 

There are LQCD data for some of the sixth order fluctuations as well. We have also computed the sixth order diagonal 
susceptibilities in PQM shown in Fig 5. The highest power of the chemical potential that appears in the SB formula 
for pressure is four (see (16)). Thus in the high temperature SB limit the sixth order susceptibilities go to zero. In 
general, these susceptibilities are very sensitive to the transition region, particularly Xe an( i Xe which stay non-zero 
only in the temperature range from 0.8T C to 1.1T C . Unlike the case of the fourth order diagonal susceptibilities where 
all of them exhibited a peak in the crossover region, the sixth order ones oscillate between positive and negative values 
passing through zero once in the transition region. Such oscillatory behaviour has also been noted in the case of 
PNJL studies [34,36]. We have also plotted the normalised susceptibilities by taking ratios with respect to the second 
order ones. As expected, these ratios approach unity in the low T limit both in PQM and HRGM. While x& IX2 
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Figure 5: On the left are plots of sixth order susceptibilities and on the right are plots of these susceptibilities 
normalized by their second order susceptibilities in ModelHotQCD (in red) and ModelWB (in blue). The lattice data 
in black are obtained using p4 action with N T = 4 by the HotQCD collaboration [31]. In all the curves the high T SB 
limit is indicated by a thick black line while the dashed black line at low T shows the HRGM value. 
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stays at unity upto almost 0.8T C , Xq /xi ancl xf/xf show deviations from unity earlier. Once again the PQM model 
predictions for the ratios are found to interpolate well between the HRGM results for low T and those of an ideal gas 
of three massless quarks in the high temperature limit. 

At zero chemical potential in the chiral limit, the singular behaviour of the quadratic and higher order baryon 
number fluctuations are supposed to be controlled by the O (4) symmetry group with the scaling behaviour [49] 



T -T r 



(25) 



where a ~ —0.25. Thus with higher order, the singular property of the baryon number fluctuations in the vicinity of 
T c grow. Similar behaviour is also expected for elctric charge fluctuations [50]. Our study with PQM for (2 + 1) flavors 
with physical light quark masses show that the behaviour of the baryon number fluctuations in the transition region is 
still controlled by (25) to a certain extent as also seen in LQCD. While the quadratic fluctuations grow monotonically 
across the transition region, the quartic fluctuations seem to show a cusp like behaviour in the crossover regime. The 
singular nature grows even stronger in the case of the sixth order fluctuations and we find oscillatory behaviour. 



4 Summary and Conclusion 

We shall now summarise our present study. In our earlier work [41], we had studied the influence of the frequently 
omitted fcrmionic vacuum fluctuations on the thermodynamics of the (2 + 1) flavor PQM and showed that the inclusion 
of this vacuum term leads to considerable improvement of the comparison between the model and LQCD. In this 
study we continue our investigation of the vacuum term in the PQM model and focus particularly on the generalised 
susceptibilities of the conserved charges B, Q and S upto order 6 for two different parameter sets of the model as 
obtained earlier [41], ModelHotQCD and ModelWB and compare our results with those of LQCD from HotQCD and 
WB groups. In this study we restrict our investigation to zero chemical potentials where LQCD data is currently 
available. 

The second order diagonal susceptibilities for B and S of ModelWB are in excellent agreement with WB lattice 
data, whereas those of ModelHotQCD fall short by ~ 7% as compared to the HotQCD data in the high T regime. In 
case of X2 > although in the high T limit there is good agreement with LQCD, in the low T limit, the ModelWB falls 
considerably short to that of the WB continuum estimate. This can be attributed to the fact that pionic fluctuations 
contribute dominantly to charge fluctuations in the low T limit which can not be captured presently in this mean field 
level analysis. The shape of the second order correlations nicely agree with that of LQCD: the BQ correlator shows 
a peak, although ModelHotQCD shows an enhanced peak by about 40% as compared to that of ModelWB. On the 
other hand, both Xii an d — Xif rises monotonically across the transition regime, with the model predictions falling 
short by ~ 30% to the HotQCD data in the transition region. 

The behaviour of the baryon number fluctuations in the transition regime are controlled by the corresponding- 
singular behaviour (25) in the chiral limit. While the second order diagonal susceptibilities show a monotonic behaviour, 
the fourth order ones tend to peak around T <~ T c and the sixth order diagonal susceptibilities oscillate once from 
positive to negative values across the crossover temeperature. In [50] it was argued that the electric charge fluctuations 
have similar features as that of the baryon number fluctuations. We confirm this expectation for the (2 + 1) flavor 
PQM in this work. In case of the fourth order fluctuations, the models do not show as pronounced a peak as lattice 
data, especially for x%-< where the slow melting of the strange condensate slows down the transition in the strange sector 
and produces a broad peak. For the sixth order fluctuations, the amplitude of the oscillation in the transition region is 
smaller compared to that of LQCD. The reasons behind these mismatch in case of fourth and sixth order fluctuations 
could be twofold: Firstly, currently LQCD data available for these quantities are only on a coarse lattice with a heavy 
pion. One has to wait for the continuum measurement with physical pion mass to make definite statements. Secondly, 
fluctuations of the field configurations could be playing a more significant role for these higher order susceptibilities 
than in case of the quadratic ones. Model computations beyond the mean field should be undertaken to verify these 
effects. We have also presented the model predictions for all the possible fourth order correlations for which there is 
no LQCD data yet. We find a sharp spike in case of the BQ correlators in the transition region while for BS and 
QS correlators, there is a sharp rise around the crossover followed by a gradual approach to the SB limit that is a 
characteristic of the slow melting of the strange condensate. There is also a hint of a double peak in case of ModelWB 
for these quantities. All the three correlators of the kind ^ B Q S show a similar behaviour, there is a sharp rise around 
T c followed by a gradual melting. 
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We have also presented results for the ratios of susceptibilities. In the low T limit, we showed that because of a 
nice factorisation in the expression of pressure in HRGM and PQM [48], the ratios become insensitive to the mass 
spectrum and thus as long as the quantum numbers of the relevant lightest degree of freedom match, one expects 
PQM to agree with HRGM. Also, on taking the ratio unknown normalisation effects cancel in case of LQCD data. 
Thus overall, the ratios of susceptibilities are in a way better suited for comparison between lattice and model results 
than the susceptibilities themselves. We find that the model prediction for the ratio of the cumulants are in far better 
agreement with lattice data. Both PQM as well as LQCD show a monotonic behaviour for X\? IX2 an d — Xif /xf 
while — XnVxf an d xfi /X2 exhibit a sharp change in the value in the crossover region. PQM also produces a sharp 
peak and a plateau for the ratios Xn /X2 an d Xn /X2 respectively as seen in LQCD. All the fourth and sixth order 
normalised diagonal susceptibilities start from unity at low T as also in HRGM and attain the SB limit in the high T 
limit. In the transition region, while the ratios of fourth to second order show a peak, the normalised sixth order ones 
show oscillatory behaviour. 

We conclude our investigation of the susceptibilities of the PQM model at zero chemical potentials and their 
comparison with available LQCD data. Overall it is encouraging to note the good agreement between the PQM model 
and LQCD data. The model successfully reproduces many of the features that are seen in LQCD data in the transition 
regime. Currently, continuum estimates for xi > X?; xf an d — xfiVxf physical pion masses are available from 
WB group. There is very good agreement between PQM and these quanitites except for xi where in the low T limit, 
the model prediction is much suppressed compared to LQCD because in the current computation at the mean field 
level pionic fluctuations that contribute dominantly to this quantity have not been included. In case of the correlations 
and higher order fluctuations, LQCD and model match qualitatively although the agreement is not as good as the 
earlier ones. Since fluctuations are expected to play an important role in these cases, model computations beyond 
the mean field might improve the agreement with LQCD. Also one has to wait for LQCD data from finer lattices 
to have a clear picture. We also found that the ratios of susceptibilities show better agreement with LQCD data 
than the susceptibilities themselves as unknown normalisation effects in LQCD data get cancelled. Also, their low 
T values approach HRGM as dependence on the mass spectrum gets cancelled. Thus, our study of these generalised 
susceptibilities and their overall good agreement with lattice data provides a strong basis for the use of PQM as an 
effective model in the study of the QCD phase diagram. In the future, we would like to compute these generalised 
susceptibilities on the entire (T — [i) plane (similar studies in the baryon sector alone have been already performed 
in [40]) that will provide us valuable insights regarding the choice of the susceptibilities best suited to search for the 
CEP. We also plan to study the effect of including fluctuations on the susceptibilities. 
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